Neural Spike Train Synchronisation Indices: Definitions, Interpretations and Applications.

A comparison of previously defined spike train syncrhonization indices is undertaken within a stochastic point process framework. The second order cumulant density (covariance density) is shown to be common to all the indices. Simulation studies were used to investigate the sampling variability of a single index based on the second order cumulant. The simulations used a paired motoneurone model and a paired regular spiking cortical neurone model. The sampling variability of spike trains generated under identical conditions from the paired motoneurone model varied from 50% { 160% of the estimated value. On theoretical grounds, and on the basis of simulated data a rate dependence is present in all synchronization indices. The application of coherence and pooled coherence estimates to the issue of synchronization indices is considered. This alternative frequency domain approach allows an arbitrary number of spike train pairs to be evaluated for statistically significant differences, and combined into a single population measure. The pooled coherence framework allows pooled time domain measures to be derived, application of this to the simulated data is illustrated. Data from the cortical neurone model is generated over a wide range of firing rates (1 - 250 spikes/sec). The pooled coherence framework correctly characterizes the sampling variability as not significant over this wide operating range. The broader applicability of this approach to multi electrode array data is briefly discussed.

Evaluation of Different Signal Processing Methods in Time and Frequency Domain for Brain-Computer Interface Applications.

Brain-computer interface (BCI) has been widely introduced in many medical applications. One of the main challenges in BCI is to run the signal processing algorithms in real-time which is challenging and usually comes with high processing unit costs. BCIs based on motor imagery task are introduced for severe neurological diseases especially locked-in patients. A common concept is to detect one's movement intention and use it to control external devices such as wheelchair or rehabilitation devices. In real-time BCI, running the signal processing algorithms might not always be possible due to the complexity of the algorithms. Moreover, the speed of the affordable computational units is not usually enough for those applications. This study evaluated a range of feature extraction methods which are commonly used for such real-time BCI applications. Electroencephalogram (EEG) and Electrooculogram (EOG) data available through IEEE Brain Initiative repository was used to investigate the performance of different feature extraction methods including template matching, statistical moments, selective bandpower, and fast Fourier transform (FFT) power spectrum. The support vector machine (SVM) was used for classification. The result indicates that there is not a significant difference when utilizing different feature extraction methods in terms of movement prediction although there is a vast difference in the computational time needed to extract these features. The results suggest that computational time could be considered as the primary parameter when choosing the feature extraction methods as there is no significant difference between the results when different features extraction methods are used.

Muscle co-activity tuning in Parkinsonian hand movement: disease-specific changes at behavioral and cerebral level.

We investigated simple directional hand movements based on different degrees of muscle co-activity, at behavioral and cerebral level in healthy subjects and Parkinson's disease (PD) patients. We compared "singular" movements, dominated by the activity of one agonist muscle, to "composite" movements, requiring conjoint activity of multiple muscles, in a center-out (right hand) step-tracking task. Behavioral parameters were obtained by EMG and kinematic recordings. fMRI was used to investigate differences in underlying brain activations between PD patients (N = 12) and healthy (age-matched) subjects (N = 18). In healthy subjects, composite movements recruited the striatum and cortical areas comprising bilaterally the supplementary motor area and premotor cortex, contralateral medial prefrontal cortex, primary motor cortex, primary visual cortex, and ipsilateral superior parietal cortex. Contrarily, the ipsilateral cerebellum was more involved in singular movements. This striking dichotomy between striatal and cortical recruitment vs. cerebellar involvement was considered to reflect the complementary roles of these areas in motor control, in which the basal ganglia are involved in movement selection and the cerebellum in movement optimization. Compared to healthy subjects, PD patients showed decreased activation of the striatum and cortical areas in composite movement, while performing worse at behavioral level. This implies that PD patients are especially impaired on tasks requiring highly tuned muscle co-activity. Singular movement, on the other hand, was characterized by a combination of increased activation of the ipsilateral parietal cortex and left cerebellum. As singular movement performance was only slightly compromised, we interpret this as a reflection of increased visuospatial processing, possibly as a compensational mechanism.

The influence of central neuropathic pain in paraplegic patients on performance of a motor imagery based Brain Computer Interface.

OBJECTIVE: The aim of this study was to test how the presence of central neuropathic pain (CNP) influences the performance of a motor imagery based Brain Computer Interface (BCI). METHODS: In this electroencephalography (EEG) based study, we tested BCI classification accuracy and analysed event related desynchronisation (ERD) in 3 groups of volunteers during imagined movements of their arms and legs. The groups comprised of nine able-bodied people, ten paraplegic patients with CNP (lower abdomen and legs) and nine paraplegic patients without CNP. We tested two types of classifiers: a 3 channel bipolar montage and classifiers based on common spatial patterns (CSPs), with varying number of channels and CSPs. RESULTS: Paraplegic patients with CNP achieved higher classification accuracy and had stronger ERD than paraplegic patients with no pain for all classifier configurations. Highest 2-class classification accuracy was achieved for CSP classifier covering wider cortical area: 82±7% for patients with CNP, 82±4% for able-bodied and 78±5% for patients with no pain. CONCLUSION: Presence of CNP improves BCI classification accuracy due to stronger and more distinct ERD. SIGNIFICANCE: Results of the study show that CNP is an important confounding factor influencing the performance of motor imagery based BCI based on ERD.

Usefulness of intermuscular coherence and cumulant analysis in the diagnosis of postural tremor.

OBJECTIVE: To investigate the potential value of two advanced EMG measures as additional diagnostic measures in the polymyographic assessment of postural upper-limb tremor. METHODS: We investigated coherence as a measure of dependency between two EMG signals, and cumulant analysis to reveal patterns of synchronicity in EMG activity in muscle pairs. Eighty datasets were analyzed retrospectively, obtained from four groups: essential tremor (ET), Parkinson's disease (PD), enhanced physiological tremor (EPT), and functional tremor (FT). RESULTS: Intermuscular coherence was highest in the PD group (0.58), intermediate in FT (0.43) and ET (0.40), and weakest in EPT (0.16) (p=0.002). EPT patients could be distinguished by low coherence: coherence <0.18 in the wrist+elbow extensors differentiates EPT in this sample with a sensitivity of 86% and specificity of 84%. Cumulant analysis showed predominantly alternating activity between wrist and elbow extensor in ET patients, while a more synchronous pattern was predominant in PD, EPT and FT (p=0.008). EMG activity in wrist and elbow flexors tended to be more synchronous in PD (p=0.059). CONCLUSION: Our results suggest that coherence and cumulant analysis may be of additional value in the diagnostic work-up of postural tremor. SIGNIFICANCE: These additional measures may be helpful in diagnosing difficult tremor cases.

The motor cortex drives the muscles during walking in human subjects.

Indirect evidence that the motor cortex and the corticospinal tract contribute to the control of walking in human subjects has been provided in previous studies. In the present study we used coherence analysis of the coupling between EEG and EMG from active leg muscles during human walking to address if activity arising in the motor cortex contributes to the muscle activity during gait. Nine healthy human subjects walked on a treadmill at a speed of 3.5­4 km h(-1). Seven of the subjects in addition walked at a speed of 1 km h(-1). Significant coupling between EEG recordings over the leg motor area and EMG from the anterior tibial muscle was found in the frequency band 24­40 Hz prior to heel strike during the swing phase of walking. This signifies that rhythmic cortical activity in the 24­40 Hz frequency band is transmitted via the corticospinal tract to the active muscles during walking. These findings demonstrate that the motor cortex and corticospinal tract contribute directly to the muscle activity observed in steady-state treadmill walking.

A portable gait assessment tool to record temporal gait parameters in SCI.

BACKGROUND: The ability to objectively analyze gait in a clinical environment is challenging due to time, space and cost constraints. This study investigated the validity of a portable gait assessment tool in objectively assessing the temporal gait parameters in subjects with spinal cord injury. The portable gait assessment tool consisted of a pair of insoles instrumented with force sensing resistors that were strategically positioned over the sole of each foot. AIM: To demonstrate the validity of the gait assessment tool by assessing the change in walking ability in incomplete spinal cord injured (ISCI) subjects, who participated in a robot-assisted gait training program. METHODS: Eighteen subjects with either an acute or chronic ISCI participated in this study (age range 26-63 years). Each subject participated in a robot assisted gait training programme for 6 weeks. Assessments were performed using the gait assessment tool before during and after the intervention. RESULTS: The gait assessment tool showed greater sensitivity to the change in the subject's gait, when compared to clinical assessments such as the walking index in spinal cord injury (WISCI II). Subjects with an acute ISCI showed a statistically significant (p<0.05) change in temporal gait parameters within the first 3 weeks of training. DISCUSSION AND CONCLUSION: This study for the first time has used the gait assessment tool in an ISCI population and has demonstrated that gait parameters can be measured and changes can be quantified within a clinical environment. The statistically significant changes during the first 3 weeks of training may indicate that an effective dose of robotic training can be administered within a relatively short period in ISCI subjects during the acute phase.

Development of quantitative and sensitive assessments of physiological and functional outcome during recovery from spinal cord injury: a clinical initiative.

The ability to detect physiological changes associated with treatments to effect axonal regeneration, or novel rehabilitation strategies, for spinal cord injury will be challenging using the widely employed American Spinal Injuries Association (ASIA) impairment scales (AIS) for sensory and motor function. Despite many revisions to the AIS standard neurological assessment, there remains a perceived need for more sensitive, quantitative and objective outcome measures. The purpose of Stage 1 of the Clinical Initiative was to develop these tools and then, in Stage 2 to test them for reliability against natural recovery and treatments expected to produce functional improvements in those with complete or incomplete spinal cord injury (SCI). Here we review aspects of the progress made by four teams involved in Stage 2. The strategies employed by the individual teams were (1) application of repetitive transcranial magnetic stimulation (rTMS) to the motor cortex in stable (chronic) SCI with intent to induce functional improvement of upper limb function, (2) a tele-rehabilitation approach using functional electrical stimulation to provide hand opening and grip allowing incomplete SCI subjects to deploy an instrumented manipulandum for hand and arm exercises and to play computer games, (3) weight-assisted treadmill walking therapy (WAT) comparing outcomes in acute and chronic groups of incomplete SCI patients receiving robotic assisted treadmill therapy, and (4) longitudinal monitoring of the natural progress of recovery in incomplete SCI subjects using motor tests for the lower extremity to investigate strength and coordination.

Kinematic modelling of a robotic gait device for early rehabilitation of walking.

Rehabilitation of walking is an essential element in the treatment of incomplete spinal cord injured (SCI) patients. During the early post injury period, patients find it challenging to practice upright walking. Simulating stepping movements in a supine posture may be easier and promote earlier rehabilitation. A robotic orthotic device for early intervention in spinal cord injury that does not require the patient to be in an upright posture has been modelled. The model comprises a two-bar mechanical system that is configured and powered to provide limb kinematics that approximate normal overground walking. The modelling work has been based on gait analysis performed on healthy subjects walking at 50 per cent, 75 per cent, and 100 per cent of normal cadence. Simulated angles of hip, knee, and ankle joints show a comparable range of motion (ROM) to the experimental walking data measured in healthy subjects. The model provides operating parameters for a prospective recumbent gait orthosis that could be used in early walking rehabilitation of incomplete SCI patients.

Human behavior integration improves classification rates in real-time BCI.

Brain-computer interfaces (BCI) offer potential for individuals with a variety of motor and sensory disabilities to interact with their environment, communicate and control mobility aids. Two key factors which affect the performance of a BCI and its usability are the feedback given to the participant and the subject's motivation. This paper presents the results from a study investigating the effects of feedback and motivation on the performance of the Strathclyde Brain Computer Interface. The paper discusses how the performance of the system can be improved by behavior integration and human-in-the-loop design.

Enhancement of a virtual reality wheelchair simulator to include qualitative and quantitative performance metrics.

The increasing importance of inclusive design and in particular accessibility guidelines established in the U.K. 1996 Disability Discrimination Act (DDA) has been a prime motivation for the work on wheelchair access, a subset of the DDA guidelines, described in this article. The development of these guidelines mirrors the long-standing provisions developed in the U.S. In order to raise awareness of these guidelines and in particular to give architects, building designers, and users a physical sensation of how a planned development could be experienced, a wheelchair virtual reality system was developed. This compares with conventional methods of measuring against drawings and comparing dimensions against building regulations, established in the U.K. under British standards. Features of this approach include the marriage of an electromechanical force-feedback system with high-quality immersive graphics as well as the potential ability to generate a physiological rating of buildings that do not yet exist. The provision of this sense of "feel" augments immersion within the virtual reality environment and also provides the basis from which both qualitative and quantitative measures of a building's access performance can be gained.

A periodogram-based test for weak stationarity and consistency between sections in time series.

In one approach to spectral estimation, a sample record is broken into a number of disjoint sections, or data is collected over a number of discrete trials. Spectral parameters are formed by averaging periodograms across these discrete sections or trials. A key assumption in this approach is that of weak stationarity. This paper describes a simple test that checks if periodogram ordinates are consistent across sections as a means of assessing weak stationarity. The test is called the Periodogram Coefficient of Variation (PCOV) test, and is a frequency domain test based on a technique of spectral analysis. Application of the test is illustrated to both simulated and experimental data (EMG, physiological tremor, EEG). An additional role for the test as a useful tool in exploratory analysis of time series is highlighted.

Reduction of common motoneuronal drive on the affected side during walking in hemiplegic stroke patients.

OBJECTIVE: The objective of this study was to use motor unit coupling in the time and frequency domains to obtain evidence of changes in motoneuronal drive during walking in subjects with stroke. METHODS: Paired tibialis anterior (TA) EMG activity was sampled during the swing phase of treadmill walking in eight subjects with unilateral stroke. RESULTS: On the unaffected side, short-term synchronization was evident from the presence of a narrow central peak in cumulant densities and from the presence of significant coherence between these signals in the 10-25 Hz band. Such indicators of short-term synchrony were either absent or very small on the affected side. Instead, pronounced 10 Hz coupling was observed. CONCLUSIONS: It is suggested that reduced corticospinal drive to the spinal motoneurones is responsible for the reduced short-term synchrony and coherence in the 10-25 Hz frequency band on the affected side in hemiplegic patients during walking. SIGNIFICANCE: This is of importance for understanding the mechanisms responsible for reduced gait ability and development of new strategies for gait restoration.

Neural spike train synchronization indices: definitions, interpretations, and applications.

A comparison of previously defined spike train syncrhonization indices is undertaken within a stochastic point process framework. The second-order cumulant density (covariance density) is shown to be common to all the indices. Simulation studies were used to investigate the sampling variability of a single index based on the second-order cumulant. The simulations used a paired motoneurone model and a paired regular spiking cortical neurone model. The sampling variability of spike trains generated under identical conditions from the paired motoneurone model varied from 50% to 160% of the estimated value. On theoretical grounds, and on the basis of simulated data a rate dependence is present in all synchronization indices. The application of coherence and pooled coherence estimates to the issue of synchronization indices is considered. This alternative frequency domain approach allows an arbitrary number of spike train pairs to be evaluated for statistically significant differences, and combined into a single population measure. The pooled coherence framework allows pooled time domain measures to be derived, application of this to the simulated data is illustrated. Data from the cortical neurone model is generated over a wide range of firing rates (1-250 spikes/s). The pooled coherence framework correctly characterizes the sampling variability as not significant over this wide operating range. The broader applicability of this approach to multielectrode array data is briefly discussed.

Organization of common synaptic drive to motoneurones during fictive locomotion in the spinal cat.

The basic locomotor rhythm in the cat is generated by a neuronal network in the spinal cord. The exact organization of this network and its drive to the spinal motoneurones is unknown. The purpose of the present study was to use time (cumulant density) and frequency domain (coherence) analysis to examine the organization of the last order drive to motoneurones during fictive locomotion (evoked by application of nialamide and dihydroxyphenylalanine (DOPA)) in the spinal cat. In all cats, narrow central synchronization peaks (half-width < 3 ms) were observed in cumulants estimated between electroneurograms (ENGs) of close synergists, but not between nerves belonging to muscles acting on different joints or to antagonistic muscles. Coherence was not observed at frequencies above 100 Hz and was mainly observed between synergists. Intracellular recording was obtained from a population of 70 lumbar motoneurones. Significant short-term synchronization was observed between the individual intracellular recordings and the ENGs recorded from nerves of the same pool and of close synergists. Recordings from 34 pairs of motoneurones (10 pairs belonged to the same motor pool, 11 pairs to close synergists and 13 pairs to antagonistic pools) failed to reveal any short-lasting synchronization. These data demonstrate that short-term synchronization during fictive locomotion is relatively weak and is restricted to close synergists. In addition, coherence analysis failed to identify any specific rhythmic component in the locomotor drive that could be associated with this synchronization. These results resemble findings obtained during human treadmill walking and imply that the spinal interneurones participating in the generation of the locomotor rhythm are themselves weakly synchronized. The restricted synchronization within the locomotor drive to motoneuronal pools may be a feature of the locomotor generating networks that is related to the ability of these networks to produce highly adaptive patterns of muscle activity during locomotion.

Effects of electrically induced muscle contraction on flexion reflex in human spinal cord injury.

STUDY DESIGN: Flexion reflex study in motor complete human spinal cord injury (SCI). OBJECTIVES: To examine changes in the magnitude of the flexion reflex following functional electrical stimulation (FES) of the rectus femoris (RF) muscle. SETTING: Bioengineering Unit, University of Strathclyde, Glasgow, Scotland, UK. METHODS: The flexion reflex was evoked by electrical stimulation of the sural nerve, and was recorded in the tibialis anterior (TA) muscle. RF muscle conditioning stimulation was performed at 0.7, 1, and 2 times motor threshold ( x MT) over a range of conditioning test intervals. RESULTS: The incidence of the early component of the flexion reflex (<100 ms) was low, suggesting that this reflex component might be suppressed in SCI. The long latency flexion reflex component (>120 ms) was observed in all subjects during control conditions and following sensorimotor conditioning. FES applied to the RF muscle (above and below MT) in the main induced a significant early and long lasting depression of the long latency flexion reflex. CONCLUSION: The depression of the flexion reflex was a result of multisensory actions on flexion reflex pathways resulting from the direct and indirect (mechanical) consequences of electrically induced muscle contraction on cutaneous and muscle afferents. Our findings emphasize the importance of sensory feedback mechanisms in modulating flexion reflex excitability, and highlight the need for rehabilitation professionals to consider the central actions of FES-induced afferent feedback when incorporating FES into a rehabilitation program. SPONSORSHIP: State Scholarships Foundation (IKY) of Hellas.

Reduction of common synaptic drive to ankle dorsiflexor motoneurons during walking in patients with spinal cord lesion.

It is possible to obtain information about the synaptic drive to motoneurons during walking by analyzing motor-unit coupling in the time and frequency domains. The purpose of the present study was to compare motor-unit coupling during walking in healthy subjects and patients with incomplete spinal cord lesion to obtain evidence of differences in the motoneuronal drive that result from the lesion. Such information is of importance for development of new strategies for gait restoration. Twenty patients with incomplete spinal cord lesion (SCL) participated in the study. Control experiments were performed in 11 healthy subjects. In all healthy subjects, short-term synchronization was evident in the discharge of tibialis anterior (TA) motor units during the swing phase of treadmill walking. This was identified from the presence of a narrow central peak in cumulant densities constructed from paired EMG recordings and from the presence of significant coherence between these signals in the 10- to 20-Hz band. Such indicators of short-term synchrony were either absent or very small in the patient group. The relationship between the amount of short-term synchrony and the magnitude of the 10- to 20-Hz coherence in the patients is discussed in relation to gait ability. It is suggested that supraspinal drive to the spinal cord is responsible for short-term synchrony and coherence in the 10- to 20-Hz frequency band during walking in healthy subjects. Absence or reduction of these features may serve as physiological markers of impaired supraspinal control of gait in SCL patients. Such markers could have diagnostic and prognostic value in relation to the recovery of locomotion in patients with central motor lesions.

Classification of Wrist Movements using EEG-based Wavelets Features.

Our aim is to assess and evaluate signal processing and classification methods for extracting features from EEG signals that are useful in developing brain-computer interfaces. In this paper, we report on results of developing a method to classify wrist movements using EEG signals recorded from a subject whilst controlling a joystick and moving it in different directions. Such method could be potentially useful in building brain-computer interfaces (BCIs) where a paralysed person could communicate with a wheelchair and steer it to the desired direction using only EEG signals. Our method is based on extracting salient spatio-temporal features from the EEG signals using continuous wavelet transform. We perform principal component analysis on these features as means to assess their usefulness for classification and to reduce the dimensionality of the problem. We use the results from the PCA as means to represent the different directions. We use a simple technique based on Euclidean distance to classify the data. The classification results show that we are able to discriminate between different directions using the selected features.

1(st) order class separability using EEG-based features for classification of wrist movements with direction selectivity.

28 channel EEG data were recorded while a subject performed wrist movements in four directions. Four feature types were extracted for each channel following optimized filtering of the signals. The potential performance of each feature and channel for use in the classification of the EEG signals was analyzed by estimating the relative class overlap using a first order histogram approach. The best feature/channel configurations contained channels both that were close and far from motor areas. While the scope and depth of the study was very limited, the results do suggest more attention should be paid to non-motor areas when investigating movement related EEG.

Functional coupling of motor units is modulated during walking in human subjects.

Time- and frequency-domain analysis of the coupling between pairs of electromyograms (EMG) recorded from leg muscles was investigated during walking in healthy human subjects. For two independent surface EMG signals from the tibialis anterior (TA) muscle, coupling estimated from coherence measurements was observed at frequencies